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Abstract: We numerically investigate the Anderson transition in an effective dimension 
d (3 < d < 11) for one particle propagation in a model random&quasiperiodic potential. 
The found critical exponents are different from the standard scaling picture. We discuss 
possible reasons for this difference. 
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The Anderson transition has been intensively investigated during last years. Accord- 
ing to the scaling theory for spinless particles all states are localized for dimension d <2, 
while a transition from localization to diffusive propagation occurs for d > 2 when a 
certain hopping parameter k crosses a critical value kcr (see reviews |2|). The one 
parameter scaling theory predicts that in the localized phase {k < kcr) the typical local- 
ization length of wavefunctions diverges at the critical point as / ~ \k — kcr\^'^ ■ Above the 
critical point the dynamics is diffusive and the DC conductivity a which is proportional to 
the diffusion rate D is assumed to approach zero at the critical point as a oc D ~ \k — kcrY. 
Scaling arguments give the following relation between the exponents: s = {d — 2)v. For 
dimension d = 2 + e the e expansion theory predicts s = v = 1 + 0(e'*)|^]. In higher 
dimensions the problem was studied by many authors ^, |6|. According to the results 
presented there i/ = l/2for(i^oo while for the delocalized phase s d/2 ^ oi a 
decreases exponentially near the transition Q. 

Numerical investigations of the exponents has been restricted to d = 3 where it was 
found s = V = 1.5ib0.1 (Q and refs. therein) in agreement with the scaling relation 
between v and s. However, the applied numerical methods were quite heavy and did not 
allow to obtain a better accuracy in the determination of the exponents or to increase 
the dimension. For 2 < d < 3 and d = 4 there are only recent results for u ^ obtained 
by transfer matrix technique. While the results there seem to be in agreement with the 
scaling theory the system size was so small that the question if the thermodynamic limit 
had been reached remains open. 

An effective way to increase the number of dimensions was proposed in HI and applied 
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for investigation of transition in d = 3 Q where it was also found s ~ 1.25 and u ~ 1.5. 
The method consists in the investigation of the well known model of quantum chaos 
namely the kicked rotator model with a frequency-modulated amplitude of kicks. The 
time-dependent Hamiltonian of the model is given by : 

H = Ho{h) + Vie,t)6i{t) (1) 

where Si{t) is the periodic delta function with period 1 between kicks, n = —id/dO and 
is a periodic angle variable. Hq determines the spectrum of unperturbed energies En 
chosen randomly distributed in (0,27r). The perturbation V depends on time in a quasi 
periodic way: 

V{e,t) = -2tan-i(2A;(cos0 + Ej=!cos(0j (2) 

with d — 1 incommensurate frequencies tOj. Here 6j are initial phases and the time is 
measured in number of kicks. The Hamiltonian (|l]) can be re-written in the extended 
phase space by letting hj = —id/dOj. After that the problem becomes periodic in time 
and the eigenvalues equation for the quasi-energy eigenfunctions can be mapped to the 
usual solid-state form [|lO|, |9| : 

TnUn + kY!r Un-r = (3) 

where the J2' indicates a sum over the nearest neighbors to n on a d-dimensional lattice 
and 

r„ = -tan(l/2(K + E-=i</'i + A)) (4) 
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Here = rijUJj and A is the quasi-energy. If (pj are randomly distributed in (0, 27r) then 
Eq. becomes equivalent to the Lloyd model at the center of the band (E = 0). The 
parameter A determines only the phase shift and it is clear that the physical characteristics 
are independent on its value. Since the mapping between (|l|) and dH) is exact, it is 
possible to study the Anderson transition in d-dimensions by investigating the dynamics 
of the 1-dimensional system (|^). This gives an effective gain of order N'^^^ in numerical 
computations if N is the system size. Finally, discussing the model, we should mention 
that the presence of disorder in the expression for is crucial. Indeed, according to the 
exact mathematical results [^] in the case of pure quasiperiodic potential when in (^) 
En = ujn all states for typical irrational frequencies are exponentially localized for any d. 
The physical meaning of this result is quite clear: the classical dynamics in this case is 
integrable and variation of unperturbed actions (levels n) is restricted by invariant curves. 
However, even if only in one direction the dispersion becomes nonlinear then the classical 
dynamics can become chaotic with diffusive spreading in all n directions. In this paper 
we investigate how this diffusion is affected by quantum effects. 

Using the above approach we studied numerically the Anderson transition for integer 
3 < c? < 11 in the model (|^) - (^). The choice of frequencies was the following: for 
d = 3 we fixed ti;i^2/2vr = 1/A,1/A^ with A = 1.3247... being the real root of the cubic 
equation — x — 1 = which gives the most irrational pair [^]; for d = 4 we added 
uj?j/2tt = l/\/2 and for d > 4 we chose all frequencies randomly in the interval (0,27r). 
The size of the basis N was between 1024 and 4096. The total number of iterations 
(kicks) was usually around 10^ but in some cases close to the critical point the evolution 
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was followed up to 10'' kicks. We used from 10 to 100 realizations of disorder to suppress 
statistical fluctuations. 

A typical example of diffusive spreading over the lattice for d = 4 is shown in Fig.l. 
Here k > kcr and the second moment of the probability distribution grows linearly with 
time. At the same time the probability distribution over levels has the gaussian shape 
(see Fig. lb). This allows to determine the diffusion constant D. Usually, we extract it 
from the probability distribution since here the statistical fluctuations are lower than for 
the value obtained from the second moment D = v? Jt. However, both methods give quite 
close values. The value of D found in this way is then averaged over different realizations 
of disorder. For k < kcr the probability distribution averaged in time reachs a stationary 
exponentially localized form and at the same time the growth of the second moment 
is saturated (Fig. 2). From the obtained stationary distribution the localization length is 
determined in the two ways. One is by the square list fitting of In I'^^p = — 2771 + b with 



the localization length being / = I/7 and b some constant [12|. Another definition is via 
the participation ratio so that 7^ = X^nlV'nl^/X^nlV'nP- After that both values 7 and 7^ 
were averaged over different realizations of disorder. The inverse participation ratio l/7i 
determines another length scale which in principle can be parametrically different from /. 

The numerical results for the dependence of D, 7 and 7, on parameter k for the effective 
dimension d = 3 are presented in Figs. 3,4. To determine the scaling near the critical 
point we used the 3 parameter fit of the type 7(j) = 70 1 A; — fecrT ^-^^ D = DQ\k — kcrl'^ . 
The results of the fitting are given in the figures captures (see also Table 1) as well as the 
parameters of test. Formally the statistical error of the exponents s, v found in this 
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way is rather small (less than 1% of the value). However, the estimate of non-statistical 
errors is quite difficult since the fitting procedure is rather sensitive to the value of kcr- 
From comparison of kcr values obtained from diffusive and localized phases it can be 
estimated on the level of 5%. The values of the exponents for D and 7^ are in the good 
agreement with the results of Q (see also Our data indicate significant difference 

for the exponents v defined via 7 and 7^ for d = 3. To demonstrate the dependence near 
the critical point we fixed the value of kcr defined from Fig. 3 and show the behaviour 
in log-log scale near kcr in Fig. 4. The two parameter fit with fixed kcr shown in Fig. 4 
gives similar values of the exponents s,i^ as in the case of Fig. 3. The linear dependence 
of In 7, In on ln|l — k/kcr\ describes quite well the variation of localization length and 
diffusion in one/two orders of magnitude. 

The case with the dimension d = A was investigated in a similar way. The results are 
presented in Figs. 5,6. They definitely show stronger deviation from the scaling relation 
between exponents s. Especially pronounced is the small value of s which remains less 
than 2. 

Inspite of this deviation from the scaling relation the behaviour at the critical point 
is close to the standard expectations [||, ||, ^, ^, ||, ^. Indeed, at kcr the conductance has 
a finite critical value g* . From another side g ~ Ed where Ec D / LP' \s the Thouless 
energy, D is the diffusion coefficient and A ^ B /L'^ is the level spacing in a block of size 
L with ~ 1 being the band width. Therefore, at kcr one has D ~ Bg* /R'^^^ where 
i? L is a typical length scale. At the same time D = B? /t so that finally R'^ ~ Bg*t. 
Our numerical data for the values of k close to kcr (Figs. 7,8) show that this relation 
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works. Formal fits give 

Rd ^ with a = 1.13(o? = 3), 1.12(d = 4) close to the expected 
value. We also analyzed the decay of the average probability to stay at the origin no: 
Po{T) = ^ dt\'i/jno{t)\'^. As can be seen in Figs 7b, 8b it is characterized by a power 
law decay. The numerically obtained values of the power are not far from l/d (Figs. 
7,8). This indicates that multifractal exponents are relatively small. This fact is also 
confirmed by the rescaling of the probability distribution = IV'nP at different moments 
of time (Fig.9). It shows that after rescaling t^^'^fn and nji}!'^ the distribution remains 
approximately stationary in time, being however not exactly exponential. 

We also investigated the delocalization transition for d = 5, 11 (Figs. 10, 11). From the 
localized side the transition was very sharp (very large 70) and it is not clear how accurate 
are the critical exponents obtained in this region even if the formal statistical error is quite 
small. Namely, the fits for 7 in the localized phase give 70 = 257(28), kcr = 0.214(1), v = 
2.32(4) with = 0.18 at d = 5 and 70 = 2115(181), A;cr = 0.107(2), 1/ = 2.55(2) with 

= 0.04 at d = 11. However, from the diffusive side the transition is going in a rather 
smooth way with the exponent s close to 2 being quite different from the expectations of 
scaling theory (see Figs. 10,11). The values of the exponents for different dimensions are 
presented in the Table 1. It definitely shows that the scaling relation s = {d — 2)u does 
not work at all. Contrary to that our numerical data indicate that for d S> 1 the exponent 
s approaches to its limiting value s ^ 2. Finally, let us mention that the dependence of 
kcr on d is quite close to kcr ^ ^/ d- This type of behaviour can be expected since in the 
kick potential all frequencies are mixed only if k Z]j=i cos{ujjt) ~ kcrd ~ 1. 

In conclusion, we studied the Anderson transition in a model random & quasiperiodic 
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potential with effective dimension d > 3. This model demonstrates quite many features 
which are the same as for the standard Anderson transition in a disordered d-dimensional 
potential. For d = 2 all states are localized and the localization length grows exponentially 
with the decrease of disorder [^] . For d = 3 the model has a transition from localization to 
diffusion with the critical exponents close to expected [^. However, for higher dimensions 
the exponents strongly deviate from the expected scaling relation s = (d — 2)v. Contrary 
to that our numerical data show that for d^ 1 one has s ~ 2. It is possible to give the 
following argument supporting s = 2. For d 3> 1 the critical value of the coupling goes to 
zero kcr ~ l/d. Therefore, the change of action is governed by the equation 

For d going to infinity this sum gives the real diffusive process in which n? ~ k'^t ~ Dt. 
Due to that it is in some sense natural to expect that asymptotically for d ^ 1 the 
exponent s approaches to 2. However, further investigations are required to conclude 
whether the behaviour of the scaling exponents is or isn't a peculiarity of the model 
under investigation. In our opinion, for the classical model chaos and diffusion in all 
directions can appear even if the motion is nonlinear only in one direction. Due to that 
we think that the quantum dynamics of the above model should be quite similar to a real 
disordered system in dimension d. 
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Laboratoire de Physique Quantique, UMR C5626 du CNRS, Univ. P.Sabatier, for the 
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(D.L.S.) thanks the University at Como for hospitality during different stages of this work. 
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Table 1: 



d 



s 



(d-2)z/(7) - s 



3 1.25(1) 1.71(6) 0.46 



2.37(1) 1.12 



4 1.52(5) 2.59(2) 3.66 



2.53(1) 3.54 



5 2.04(3) - 



2.32 



4.92 



11 1.87(1) - 



2.55 



21.0 



Figure 1: (a) Behaviour of the second moment on time for k — 0.33, d — 4: (diffusive 
side) . Dashed hne is the fitted diffusion coefficient D — 0.0205. Fourier basis is N — 
1024. (b) Logarithm of the averaged distribution function between t = 0.995 • 10*^ 
and t — 10^ (t=number of kicks). Full line is the best fit gaussian which gives a 
diffusion coefficient D — 0.0229. On x-axis we put n^. 

Figure 2: (a) Behaviour of the second moment on time for k = 0.24, d = 4 (localized 
side). Basis is N — 1024. (b) Logarithm of the averaged distribution function 
between t — 0.3995 • 10^ and t — 0.4 • 10^. On x-axis we put n. Line represents the 
fitted localized distribution with / = 15.4. 

Figure 3: Inverse localization length (left side) and diffusion rate (right side) for 
d — 3. Lines are the separate fits: 7 = joikcr — ky with 70 = 8.01(1), kcr = 0.509(1) 
and u = 2.37(1) with = 10-5 and D = Do{k - k^Y with Dq = 2.56(2), k^r = 
0.479(1) and s = 1.25(1) with = 32.6. A similar fit for 7^ = -foi{kcr - k^' gives 
7oi = 3.7(3), kcr = 0.489(4), i/^ = 1.71(6) and = 15.4. Here 10 - 100 random 
configurations have been iterated up to 10^ kicks. Data errors are within the symbol 



size. 
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Figure 4: d — 3. (a) Logarithm of diffusion rate vs ln{k/kcr — 1) where kcr is 
extracted from the separate fit (Fig. 3). Straight hne is the best fit hne with slope 
s — 1.248(5) and — 10-6. (b) Logarithm of the inverse locahzation length vs 
ln{l — k/kcr) where kcr is extracted from the separate fit (Fig.3). Straight hne 
is the best fit line with slope u = 2.374(9) and = 33.8. (c) Logarithm of the 
participation ratio 7j vs ln{l — k/kcr) where kcr — 0.489(4) is extracted in a similar 
way from a three parameter fit. Here the straight line has slope Ui = 1.71(6) and 
= 16.3 Here 10 — 100 random configurations have been iterated up to 10^ kicks. 
Errors are within the symbol size. 



Figure 5: Inverse localization length (left side) and diffusion rate (right side) for 
d = 4. Lines are the separate fits: 7 = ■jo{kcr — ky with 70 = 66.(10), kcr = 0.306(2) 
and u = 2.53(1) with = 10-2 and D = Do{k - kerf with Dq = 2.29(2), kcr = 
0.283(3) and s = 1.52(5) with = HI- A similar fit for 7^ = ^Qi{kcr - k)"' gives 
7oi = 85(28), kcr = 0.305(5), = 2.59(2) and x^ = 4.5. Here 10 - 100 random 
configurations have been iterated up to 10^ kicks. Data errors are within the symbol 
size. 
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Figure 6: d — A. (a) Logarithm of diffusion rate vs ln{k/kcr — 1) where kcr is 
extracted from the separate fit (Fig. 5). Straight hne is the best fit hne with slope 
s — 1.519(4) and — US. (b) Logarithm of the inverse locahzation length vs 
ln{l — k/kcr) where kcr is extracted from the separate fit (Fig 5). Straight line is 
the best fit line with slope v = 2.534(5) and = 91.81. (c) Logarithm of the 
participation ratio 7j vs ln{l — k/kcr) where kcr — 0.305(5) is extracted in a similar 
way from a three parameter fit. Here the straight line has slope Ui = 2.59(2) and 
= 4.5 Here 10 — 100 random configurations have been iterated up to 10^ kicks. 
Errors are within the symbol size. 

Figure 7: d — 3 Study near the critical point at k=0.48. (a) Behaviour of R''' — 
J2n\iT'\'^\'4'n\'^ in time. Best fit (dashed line) gives ~ ^i.i3(3)_ ^j^-j Behaviour 
of the logarithm of the integrated probabihty to stay Po{T) as a function of the 
logarithm of time T. Here is Po{T) = ^ Jq dt\ipno{t)f ^^'i the initial state vector 
is ilj{t = 0) = Sn,no- Best fit (dashed line) gives Po{T) - T~^-n3)_ Qne random 
configuration has been considered. Basis is N=4096. 

Figure 8: The same as Fig. 7 but for d=4 near the critical point at k=0.27. Results 
from fit gives R'^ - t^-^^Cs) p^^^"^ _ ^.-.26(3) 
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Figure 9: Rescaling of the distribution function in the critical region {k — 0.27) 
for d = 4. (a) The averaged (over 10^ kicks ) probabihty distribution over the 
unperturbed basis taken at three different times ti,2 ti,3 ti where ti — 10^, 
shows that as the time goes on the distribution increases its size, (b) The same as 
(a) but in the rescaled variables t^^'^fn and n/t^^'^ (since at the critical point |n|'^ ~ t 
with d = 4). 



Figure 10: d — (a) Diffusion rate vs /c, 3 parameters fit (full hne) gives : 
D = 2.56(6) {k - .180(2))2°4(3) ^ ^ 35 (b) logarithm of diffusion rate 

vs ln{k/kcr — 1) where kcr is extracted from (a). Best fit (full line) has the slope 
s = 2.039(9) 



Figure 11: d — 11. (a) Diffusion rate vs /c, 3 parameters fit gives : D — 1.35(1) {k — 
.0924(5))^'^^*^^^ with = 310; (b) logarithm of diffusion rate vs ln{k/kcr — 1) where 
kcr is extracted from (a) Best fit (full line) has the slope s — 1.99(5) 
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